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Abstract 

In this paper, optical and electronic properties of diamond shaped graphene quantum dots 
(DQDs) have been studied by employing large-scale electron-correlated calculations. The 
computations have been performed using the 7r-electron Pariser-Parr-Pople model Hamil¬ 
tonian, which incorporates long-range Coulomb interactions. The influence of electron- 
correlation effects on the ground and excited states has been included by means of the 
configuration-interaction approach, used at various levels. Our calculations have revealed 
that the absorption spectra are red-shifted with the increasing sizes of quantum dots. It 
has been observed that the first peak of the linear optical absorption, which represents the 
optical gap, is not the most intense peak. This result is in excellent agreement with the 
experimental data, but in stark contrast to the predictions of the tight-binding model, ac¬ 
cording to which the first peak is the most intense peak, pointing to the importance of 
electron-correlation effects. 


1 


I. INTRODUCTION 


Graphene, a two-dimensional, one atom thick layer of graphite, with carbon atoms ar¬ 
ranged in a honeycomb lattice, has attracted enormous attention among researchers in re¬ 
cent years due to the possibilities of appealing applications in the held of nanoelectronicsP 
However, a major drawback of pure graphene from the viewpoint of electronic devices, in 
general, and opto-electronic devices, in particular, is its zero band gap. This problem has 
stimulated tremendous amount of experimental and theoretical efforts in attempting to for¬ 
mulate techniques to introduce a band gap in grapheneP It has been observed that reducing 
the dimensionality of graphene opens up the band gap on account of quantum confinement. 
One-dimensional periodic graphene nanostructures such as graphene nanoribbons (GNRs) 
have band gaps ranging from zero (metallic) to rather large values (semiconducting), de¬ 
pending on the width of the ribbon, and the nature of edge termination. 3 Opening up the 
band-gap further, by reducing the dimension of GNRs, has been made possible by fabrica¬ 
tion of stable zero-dimensional graphene quantum dots (GQDs).® 6 The band-gaps of GNRs 
are ?a0.4 eV, while the band-gaps of GQDs can be tuned up to ~3 eV by decreasing their 
size! 7 This appealing feature of GQDs enhances the prospects of utilization of such mate¬ 
rials in lasers, light-emitting diodes (LEDs), solar cells, bio-imaging sensors,® and optically 
addressable qubits in quantum information science.® Since, electronic excitation determine 
the photophysics of GQDs which is vital for all these applications, it is essential to have 
a detailed understanding of their low-lying excited states which can be probed by optical 
means. 

Because of aforesaid possibilities of applications of GQDs, significant studies, experi¬ 
mental, as well as theoretical, on the electronic and optical properties of GQDs have been 
performed lately. Experimental studies on GQDs (in the size range of 5-35 nm) by Kim et 
al.‘ have revealed that while the absorption peak energies decrease with increasing size of 
the GQDs, the photoluminescence (PL) spectra exhibit a decrease in energy as the average 
size of GQDs increases up to ~17 nm, followed by an increase in the peak energy with 
increasing average size of GQDs. They accredited this abnormal behavior of PL spectra to 
the presence of edge variations associated with size-dependent shape in GQDs. However, 
single-particle spectroscopic measurements carried out by Xu et al. U) have shown that size 
differences of GQDs do not affect the peak positions and spectral line-shapes. Experiments 
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on photoexcited GQDs have also indicated that emission intensity decreases when GQDs are 
excited to singlet states (Si, S 2 , S 3 ), while it increases sharply when they are excited to S 4 , or 
higher excited states) 11 In addition, several experimental studies have shown that the optical 
band-gap is dependent on the size of the GQDs, giving rise to different excitation/emission 
spectra as well as PL spectra of different colors)- 12 13 Further, it has been observed that PL 
behavior is strongly associated with the presence of microstructures in GQDs and hence, 
are affected by edge effects as well as emission sites. 11 Thus, it is essential to have a detailed 
knowledge of the atoms which significantly contribute to the optical band-gap. 

As far as theoretical studies are concerned, Yamijala et a/. 15 have performed a detailed 
study of the structural stability, electronic, magnetic and optical properties of rectangular 
shaped graphene quantum dots as a function of their size, as well as under the application 
of electric held, using first-principles density functional theory (DFT). However, DFT based 
calculations are known to underestimate the band gap, and provide a reasonable description 
of the excited states only when they do not exhibit significant configuration mixing. Yan 
et al. 16 have employed a tight-binding (TB) model to predict the band gap of GQDs as a 
function of their size. Theoretical calculations using the TB model have also been utilized 
to study the optical properties of hexagonal 17 and triangular graphene quantum dots as a 
function of their size and type of edge) 18 However, TB method is unreliable in predicting low- 
lying excited states because it does not include electron-electron interactions. In addition, 
calculations of optical properties of large graphene quantum dots employing first-principles 
DFT have also been performed by Schumacher. 19 Recent first principles DFT calculations 
on diamond-shaped graphene nanopatches have revealed that these systems display well- 
defined magnetic states which can be selectively tuned by the application of electric held. 20 
However, to the best of our knowledge, till date there is no existing literature (experimental 
as well as theoretical) on the optical properties of diamond shaped graphene quantum dots. 

Motivated by aforementioned theoretical and experimental studies of GQDs, in this work 
we present a systematic study of the electronic structure and the optical properties of 
diamond-shaped graphene quantum dots (DQDs) which exhibit a mixture of zigzag edges 
and armchair corners, and we hope that our studies will motivate experimentalists to explore 
optical properties of these nanostructures. For this purpose, we have employed a method¬ 
ology based upon Pariser-Parr-Pople (PPP) model Hamiltonian, 21 22 which is an effective 
7T-electron model, including long-range electron-electron interactions. We have used this 
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approach in several works in our group dealing with conjugated polymers,® 29 polyaromatic 
hydrocarbons,® 0 31 graphene nanoribbons,® 33 and graphene nano disks.® PPP model has an 
advantage over the TB model in that it incorporates long-range Coulomb interactions among 
the 7r-electrons, essential for taking into account influence of electron correlation effects. Fur¬ 
ther, it considers the interactions of 7r-electrons with a minimal basis, therefore, as compared 
to ab initio approaches, it yields highly accurate results with fewer computational resources. 
In this work, we present theoretical calculations of the electronic structure and linear opti¬ 
cal absorption spectra of DQDs of varying sizes employing a configuration-interaction (Cl) 
methodology) 2 ®® so as to account for electron-correlation effects in their ground and excited 
states. As far as experiments are concerned, it is impossible to synthesize bare graphene 
quantum dots of high symmetry, because, due to the dangling bonds, edges will undergo 
significant reconstruction, leading to distorted shapes. Nevertheless, several polycyclic aro¬ 
matic hydrocarbons (PAHs) have been synthesized which are nothing but graphene quantum 
dots of high symmetry, but with edges passivated by hydrogen atoms,® a few of which we 
had studied in earlier works) 30 31 Of the quantum dots considered here, hydrogen passivated 
counterpart of DQD with 16 carbon atoms (DQD-16, henceforth) is called pyrene, while that 
of DQD with 30 carbon atoms (DQD-30) is known as dibenzo[bc,kl]coronene, both of which 
have been well-studied in the chemical literature. 35 A large number of experimental measure¬ 
ments of optical absorption of pyrene in vapor,®solution, 37 40 and matrix isolated phased®® 
have been performed, and our results on DQD-16 are in excellent agreement with them. Clar 
and Schmidt® measured the gas-phase absorption spectrum of dibenzo[bc,kl]coronene, and 
our calculations on DQD-30 are in very good agreement with the experimental results. We 
also computed the absorption spectrum of next larger quantum dot DQD-48, whose struc¬ 
tural properties have been studied theoretically by several authors,® 49 but it has not been 
synthesized as yet. 

Theoretically, Canuto et a/.,® computed the absorption spectrum of pyrene employing 
the intermediate neglect of differential overlap (INDO/S) semi-empirical quantum mechan¬ 
ical technique along with singles configuration interaction (SCI) method, while Guclipati et 
a/.,® calculated the excitation energies and oscillator strengths of pyrene using the com¬ 
plete neglect of differential overlap (CNDO/S) model, coupled with the truncated singles 
and doubles configuration interaction (SDCI) method. Parac et al.f ^ and Malloci et a/.® 
employed the time dependent density functional theory (TDDFT) technique to compute 
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the excitation energies and photoabsorption cross-sections of pyrene, respectively. Malloci 
et a/. ® also calculated the photoabsorption cross-section of dibenzo[bc,kl]coronene using 
the TDDFT technique. Several authors have studied the structural stability of the PAH 
equivalent of DQD-48 (C4 8 Hi 8 ) by employing first-principles DFT-based methodologies.®®! 
Additionally, Karki et al , 47 also studied the variation of the optical gap with increasing size 
of the PAH clusters, while Boersma et a/.,® and Pathak et al . 4 -^ calculated their infrared 
spectra. Denis et a/.,® analysed the effect of addition of azomethine ylide on the binding 
energy of C 48 H 18 . 

Based upon our calculations, we predict the variation in the behavior of linear absorption 
spectrum with increasing size of DQDs, and our results are in significant variance with the 
predictions of the TB model. We also identify the atoms which play a significant role in the 
band gap, and thus the optical spectrum of the DQDs. 

The remainder of this paper is organized as follows. In section [H] we present a brief 
overview of the theoretical methodology adopted by us. In section H we present and 


discuss the results, followed by conclusions in section IV An Appendix representing detailed 
information about many-particle wave-functions of excited states contributing to the optical 
absorption peaks is presented at the end of the paper. 


II. COMPUTATIONAL DETAILS 

The schematic diagram of the geometry of DQDs considered in this work is given in 
FigjTJ Different DQDs can be identified by the total number of carbon atoms n, and will 
be denoted as DQD-n, henceforth. In our calculations, all quantum dots are assumed to lie 
in the x — y plane, with the shorter diagonal of the DQD assumed to be along the x-axis, 
and the longer one along the r/-axis. All carbon-carbon bond lengths and bond angles have 
been chosen as 1.4 A, and 120°, respectively. The point group of DQDs is D- 2 h , with PA^ 
being the ground state. Then, as per electric-dipole selection rules, the symmetries of the 
one-photon excited states are 1 B 2u and 

These calculations have been carried out by employing the PPP model Hamiltonian,122221 
given by 
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Figure 1: Schematic diagram of diamond shaped graphene quantum dots consisting of 16 

atoms, 30 atoms and 48 atoms, respectively. 

where cj^c^) creates (annihilates) a 7r orbital of spin a , localized on the zth carbon atom 
while the total number of electrons with spin cr on atom i is indicated by n* = X )a c \a c ia■ 
The second and third terms in Eq. [l] denote the electron-electron repulsion terms, with the 
parameters U, and Vij representing the on-site, and the long-range Coulomb interactions, 
respectively. The matrix elements tij depicts one-electron hops, which in our calculations, 
have been restricted to nearest neighbours, with the value t 0 = 2.4 eV, consistent with our 
earlier calculations on conjugated polymers,® 29 and polyaromatic hydrocarbons. 30 31 

Parameterization of the Coulomb interactions is done according to the Ohno relationship^ 


V lj = U/K iJ {l + 0.6ll7Rl j ) 1/2 (2) 

where Kij represents the dielectric constant of the system which replicates screening 
effects, U as described above is the on-site electron-electron repulsion term, and Rij is the 
distance (in A) between the ith and jth carbon atoms. In the present work, we have done 
calculations adopting both “screened parameters’^- 55 with U = 8.0 eV, Kij = 2.0 (i ^ j) 
and Kij = 1.0, as also the “standard parameters” with U = 11.13 eV, and Kij = 1.0. We 
observe that our calculations employing the screened parameters, proposed by Chandross 
and MazumdarpS are in better agreement with the experimental results, as compared to 




those performed using standard parameters, consistent with the trends observed in our 
earlier works as well. 212 ' 


The first step of our calculations is to find the self-consistent solutions at the Restricted 


Hartree-Fock (RHF) level, employing the PPP Hamiltonian (c/. Eq. [I]), using a code devel¬ 
oped in our group . 34 These solutions, in which electrons occupy the lowest energy orbitals, 


comprise the HF ground state. This is followed by correlated calculations at the quadruple 
configuration interaction (QCI) level, or at the multi-reference singles-doubles configuration 
interaction (MRSDCI) level, depending upon the size of DQD. In the QCI approach, up to 
quadruple excitations from the HF ground state are considered, and, thus, it requires sig¬ 
nificant amount of computational resources. Therefore, QCI calculations can be performed 
only for small systems (in our case, for DQD-16). For larger DQDs, MRSDCI approach has 
been employed. In MRSDCI calculations, singly and doubly excited configurations from the 
reference configurations of the selected symmetry subspace are considered while generating 
the Cl matrix . 56 — Subsequently these Cl wave functions are used to compute transition 
electric dipole matrix elements between various states, required for computing the optical 
absorption spectra. From the calculated spectra, important excited states giving rise to 
various peaks are identified, and the dominant reference configurations contributing to these 
excited states are included to enhance the new reference space. This procedure is iterated 
until the desired absorption spectrum converges to an acceptable tolerance. With the in¬ 
creasing sizes of the DQDs, the number of molecular orbitals of the DQD increases, leading 
to an increase in the size of the Cl expansion. Therefore, to make calculations feasible, the 
frozen orbital approximation was adopted for DQD-48, with the lowest two occupied orbitals 
frozen, and highest two virtual orbitals deleted, so as to retain the particle-hole symmetry. 

The formula employed for the calculation of the ground state optical absorption cross- 
section cr(u}), assumes a Lorentzian line shape 



( 3 ) 


where u denotes incident radiation frequency, e denotes its polarization direction, r is the 
position operator, a is the Ene structure constant, 0 and i denote, respectively, the ground 
and the excited states, (U;o is the frequency difference between those states, and 7 is the 
absorption line-width. 


7 







Table I: Dimensions of Cl matrices for DQDs of varying sizes, for various symmetry 
manifolds. QCI method was employed for carrying out calculations on l A g and 1 7?2« 
manifolds of DQD-16, while for rest of the calculations MRSDCI method was used. 


Number of atoms in DQD 

Dimension of Cl matrix 


1 4 

A g 

1 b 2u 

1 b 3u 

16 

73857 

126279 

142992 

30 

215919 

1564554 

1359014 

48 

237030 

5442399 

4269236 


III. RESULTS AND DISCUSSION 

In this section, we present the results obtained from Cl calculations for DQDs of varying 
sizes, ranging from DQD-16 to DQD-48. In order to acquaint the reader with the precision 
of our MRSDCI or QCI calculations, the sizes of the resultant Cl matrix for different sym¬ 
metries for the DQDs considered here are given in Table [T| QCI method was employed for 
carrying out calculations on 1 A g and 1 B 2u manifolds of DQD-16, while rest of computations 
on DQD-16, DQD-30, and DQD-48, were carried out by adopting MRSDCI methodology. 
Sizes of the Cl matrix indicate that electron-correlation effects were well accounted for in 
these calculations. 


A. Charge density distribution and optical gap 

In this section we examine the evolution of the band gap, orbital energy levels, and the 
charge densities with the size of the DQD. 


1. Charge Density 

The charge density bubble-plots for the HOMO orbital obtained by employing the TB 
model and the PPP model for DQDs of varying sizes are presented in Fig. |2j The charge 
density of the LUMO orbital is same as that of the HOMO orbital because of the electron- 










hole symmetry, and hence has not been shown. The numbering scheme of atoms for the 
different DQDs considered is also presented. 

It is observed that in case of the smallest quantum dot DQD-16, the contribution from 
atomic sites 1, 3, 14 and 16 to the charge density of HOMO orbital is maximum. These 
atoms are at the projected corners of a purely zigzag edge as is evident from Fig. [2j The 
charge density contribution from the atomic sites 4, 8, 9, and 13, which are also at the edge 
of the diamond quantum dot, is less as compared to that of the atoms mentioned earlier. 
This can be attributed to the fact that these atoms give rise to an edge which exhibits both 
zigzag as well as armchair nature. With the increasing size of the quantum dot, the zigzag 
characteristic of the edge becomes more conspicuous, leading to increased contribution of the 
atoms located on the zigzag edges. This trend is obvious from the dominant contributions to 
the charge densities of the HOMO orbital by atoms 4, 8, 23 and 27 for DQD-30, and atoms 
4, 8, 41 and 45, for DQD-48. From this it is evident that the HOMO-LUMO band-gap, and 
thus the optical properties of DQDs, can be tuned if suitable functional groups are attached 
to these atoms on the zigzag edges. This is in agreement with results obtained earlier for 
the case of graphene nano-ribbons (GNRs),® as well as triangular nanographenesP^ 

Looking at the energy-level diagrams of various DQDs presented in Fig. [3j unlike the case 
of triangular nanographenes with zigzag edges, one notices the absence of zero-energy states. 
It is also evident from the figure that the increasing size of DQDs leads to the reduction of 
the band gap, consistent with the gapless nature of infinite graphene. 


2. Optical gap 

In Table [TTJ we present the energy gap between the highest-occupied molecular orbital 
(HOMO), and the lowest-unoccupied molecular orbital (LUMO), for increasing sizes of 
DQDs obtained from the TB model, as well as from the PPP model, at the RHF level. 
Because the band gap is also the optical gap in DQDs, in the same table we also present 
the results of the optical gap of these systems obtained from our correlated electron Cl 
calculations. We note the following trends in these results: (a) the gaps decrease with the 
increasing DQD size, irrespective of the Hamiltonian or the method used, (b) gaps obtained 
from the TB model are significantly smaller as compared to those obtained from other meth¬ 
ods, (c) at the HF level standard parameter gaps are much larger than those obtained using 
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Figure 2: Charge density bubble plots of HOMO orbital obtained by employing TB and 
PPP models for DQDs of varying sizes. The numbering scheme of atoms for the different 

DQDs considered is also presented. 
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Figure 3: Energy level plots for DQDs of varying sizes, calculated using the TB model. 

screened parameter, (d) optical gaps at the Cl level are significantly redshifted as compared 
to their HF values, in the standard parameter calculations. But, for the screened parameter 
calculations, these correlation-induced shifts are small, with the Cl level gaps of DQD-16 
and DQD-30 exhibiting slight redshifts, while that of DQD-48 exhibiting a small blueshift, 
and (e) at the Cl level the gaps obtained using the two sets of PPP parameters are in 
good quantitative agreement with each other, suggesting the correctness of our correlated 
electron approach. This decrease of the gap with the increasing sizes observed for DQDs 
is in agreement with previous experimental and theoretical results obtained for graphene 
quantum dots of other shapes.®^ 


B. Linear Absorption spectrum 

In this section, we first elucidate the salient features of the linear optical spectra of DQDs 
of varying sizes computed within the framework of the independent-electron TB model and 
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Table II: HOMO-LUMO band gap for increasing size of DQDs obtained from the TB 
model and the PPP model. In case of PPP model, the gap is calculated both at the HF 
and Cl level, using the standard (Std) as well as the screened (Scr) parameters. At the Cl 

level, the gap is identified with the optical gap. 


System 

HOMO-LUMO gap (eV) HOMO-LUMO gap (eV) Optical gap (eV) 

(TB model) (PPP-HF) (PPP-CI) 



Scr 

Std 

Scr 

Std 

DQD-16 

2.14 

3.91 

7.26 

3.60 

3.74 

DQD-30 

0.89 

2.11 

4.65 

2.08 

2.31 

DQD-48 

0.34 

1.10 

2.81 

1.40 

1.61 


PPP model at the HF level, which will allow us to gauge the influence of electron-correlation 
effects in the PPP model-based Cl calculations, presented thereafter. 


1. Calculations at the tight-binding and HF level 

The absorption spectrum obtained from TB model and PPP model at HF level employing 
screened parameters (Fig. [4]) exhibits the following characteristics: 

1. The absorption spectrum is red-shifted with increase in size of the DQD, in agreement 
with quantum confinement effect. This red-shift is more pronounced at the PPP-HF 
level as compared to that obtained by the TB model. In addition, the absorption 
spectrum at the PPP-HF level is blue-shifted compared to the one computed using 
the TB model because it is well known that the HF theory overestimates energy gaps, 
as was also observed in earlier works!- 59 60 

2. The pattern of the absorption spectra at the PPP-HF level is similar to that obtained 
by the TB model. However, the absolute intensities of the peaks at the PPP-HF level 
are lesser as compared to those obtained by the TB model. 

3. The first peak at the TB and PPP-HF level is always ^/-polarized, and corresponds to 
excitation of a single electron from the HOMO ( H ) orbital to the LUMO (L) orbital. 
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This peak is also the most intense peak in the calculated spectra, which is in stark con¬ 
trast with the experimental results obtained for pyrene and dibenzo[bc,kl]coronene. 36 "— 61 

4. With the increasing size of the DQDs, the intensity of the first peak increases enor¬ 
mously as compared to the other peaks in the absorption spectrum. For example, in 
DQD-30/DQD-48, the relative intensity of the peaks starting from the second one is 
much smaller as compared to DQD-16, irrespective of the Hamiltonian employed. 

5. We also note that ^-polarized peaks are degenerate, while y-polarized peaks exhibit 
no degeneracy. For example, in case of DQD-16, the second peak is x-polarized and 
is due to degenerate excitations \H — 1 —>• L) and \H —> L + 1), while the third peak 
is y-polarized, and is due to nondegenerate excitation \H — 1 —> L + 1). 


2. Correlated-electron calculations 


We discuss the general features of the optical absorption spectra obtained from the Cl 
calculations, followed by a more detailed examination of individual DQDs. 

One observes the following general trends upon examining the absorption spectra calcu¬ 
lated by the PPP-CI approach presented in Figs. [5] [7| and the quantitative information 
about various excited states detailed in Tables |I\ | [X] of the Appendix: 


1. In agreement with the TB results, the absorption spectrum is red-shifted with the 
increasing size of the DQD. 


2. The spectra obtained from screened parameters is red-shifted as compared to that 
obtained from standard parameters. Furthermore, for the case of DQD-lfP 6 ® and 
DQD-30 , 145,61 1 screened parameter results are in overall better agreement with the ex¬ 


perimental data (c/. Table III), as compared to the standard parameter ones. 


3. The first peak in the spectrum of the DQDs is always due to the absorption of a 
y-polarized photon, causing a transition from their ground state (l 1 ^), to the 1 1 B 2 U 
excited state, and denotes the optical gap. The wave function of the l 1 B 2 U state for 
all the DQDs is dominated by the \H —> L) excitation, in agreement with the results 
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Table III: Comparison of locations of experimentally measured linear absorption peaks of 
pyrene and dibenzo[bc,kl]coronene with our PPP-CI results for DQD-16 and DQD-30, 
respectively. Computed values of the excited state energies for DQD-48 are also presented. 
DF and MI, respectively, denote dipole forbidden state, and most intense peak. Theoretical 
results of other authors are also presented for comparison. All energies are in eV units. 


System 

Symmetry 

Experimental Values 

Theory 

(Other) 

This work 





Scr Std 

DQD-16 

B 3u (DF) 

3.34, 1 33 3 .33, S3 

3.33,22 3.83,22 3.7^3 

2.82 3.09 


B211 

(optical gap) 

3.70,® 3.69, 1221221 3 . 71 ® 3.80,SU 

3.75,® 3.79,® 3.8 # 1 

3.35,22 3 . 53 EE 3 . 93) Sa 3 . 

3.60 3.74 


b 3u 

4.55,™ 40 4.62,22 

4.67, 41 44 4.60,22 

4.10,® 4.70,® 5.2d® 

4.37 4.94 


B211 

5.15p2 5 . 12 , 1221221 5.17,22 5.29,22 
5.34, 41 441 5.35,22 5.2^32 

5.00,22 5.36 ,® 5.6cJ® 

5.37 5.44 


B211 

6.32,22 6.07® 

5.80,® 5.96,® 6 . 6 d® 

6.22 6.44 


B 3u 

6.4#2 

6.04,® 6.35 ,22 6.9d® 

6.38 6.96 


u 

7.0^32 

6.45,® 7.39,® 7.5l22 

6 . 88 , 6.87 7.29, 7.33 


B 2u ( MI) 

5.15,23 5.12, 1 221221 5.17,22 5.29,22 
5.34, 4 f2il 5.35,22 5.2^33 

6.35,® 5.36 ,® 5.6CP 

5.37 6.44 

DQD-30 

B2u 

(optical gap) 

2.5^33 

2 id 35153 

2.08 2.31 


B 3u (DF) 



2.25 2.43 


B 3u 

3.61® 

3.5# r>:! 

3.45 3.83 


B. 3 „(MI) 

6 . 2 cP 

5.9# 53 

6.28 5.10 

DQD-48 

B 3u (DF) 



1.38 1.52 


B211 

(optical gap) 



1.40 1.61 


B 2 u (MI) 



2.19 5.54 
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Figure 4: Optical absorption spectrum of DQD-16, DQD-30, and DQD-48, calculated 
using the tight binding model and PPP model at HF level employing screened parameters. 

The spectrum has been broadened with a uniform line width of 0.1 eV. 

of the TB model. However, a quantitative analysis of the optical gap indicates that its 
value obtained from the TB model is much less compared to the value obtained from 
the PPP-CI approach. For the cases of DQD-lfP^® and DQD-30, 45 61 for which the 
experimental results are available, again PPP-CI value of the gap is in much better 
agreement with experimental results, than the TB model value. Therefore, we hope 
that similar experiments can be performed on DQD-48 in the future, so that our 
predicted PPP-CI values of optical gaps can be tested. 

4. The intensity of the first peak is lesser as compared to the intensity of other higher 
energy peak(s) in the spectrum, in contrast to the predictions of the TB model, and in 
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agreement with the experimental results for pyrene®'® and dibenzo[bc,kl]coroneneJ 45 61 

5. The calculated position of the first peak in the PPP-CI absorption spectrum is weakly 
dependent on the choice of the Coulomb parameters in the PPP model (standard or 
screened). However, higher energy peaks, and the character of the many-particle wave 
functions contributing to them, do depend significantly upon the choice of Coulomb 
parameters. In particular, the position of the most intense peak is drastically depen¬ 
dent upon the choice of the Coulomb parameters in the PPP-CI calculations. Thus, 
we conclude that the position of the most intense peak is strongly dependent on the 
strength of the Coulomb interactions in theses systems. 


6. The optical transition to the first excited state of B 3u symmetry is dipole forbidden 
within the PPP model, on account of the particle-hole symmetry due to the use of the 
nearest-neighbour hopping approximation. However, this symmetry is approximate in 
the real systems, and hence the transition to this state is experimentally observed as 
a weak peak. Our PPP-CI calculations predict this state to lie below the optical gap 


for DQD-16 and DQD-48, but above it for DQD-30. (cf. Tables IV V VI VII VIII 


and 0. Our prediction is in agreement with the experiments for the case of DQD-16 
when compared with the data for pyrene (cf. Table III ), however, no experimental 
results for this state are available for the larger DQDs. 


7. Wave functions of higher energy states derive significant contributions from double 
and higher level excitations, signaling the importance of electron correlation effects. 


DQD-16 


Figure [5] presents the computed linear absorption spectrum for DQD-16, obtained by em¬ 
ploying standard as well as screened parameters, while Tables IV and [V] of the Appendix, 
present the detailed quantitative data corresponding to various peaks in the computed spec¬ 
tra, and the excited states contributing to them. Our theoretical results have been compared 
with the experimental optical absorption data of pyrene (CieHio), which is nothing but hy¬ 
drogen saturated DQD-16. Because we have employed PPP model parameters used to 
describe the optical properties of aromatic hydrocarbons, the comparison between DQD-16 
and pyrene is most appropriate. A number of theoretical 50 52 and experimental studies 36 ® 


16 
















of optical absorption in pyrene have been carried out in the past, and our calculated excited 
state energies are found to to be in very good agreement with the results obtained earlier 


(cf. Table 111 


The first peak in the experimentally obtained absorption spectrum of pyrene is a weak 
one located around 3.34 eV, 37 and corresponds to the dipole forbidden B 3u state in our 
calculations, mentioned earlier. Our standard parameter value of 3.09 eV for the excitation 
energy of this state is in good agreement with the experimental value, while the screened 
parameter value of 2.82 eV underestimates it. The location of the second peak in the 
experimental spectrum (3.69-3.83 eV), which also defines the optical gap, is in excellent 
agreement both with our standard and screened parameter PPP-CI values of the optical gap, 
computed at 3.74 eV, and 3.60 eV, respectively. This peak is //-polarized and corresponds 
to lBiu state. The optical transition to the fourth excited state gives rise to the most 
intense peak experimentally observed to be in the range 5.15-5.35 eV. This result is in 
excellent agreement with our PPP-CI value 5.37 eV, obtained using the screened parameters. 


As a matter of fact, it is obvious from Table III that the agreement between the dipole- 
allowed states obtained from our screened-parameter based PPP-CI calculations, and the 
experimental measurements of Becker et al. 37 and Gudipati et al. 42 is quite remarkable both 
for peak locations, and the symmetry assignments, all the way up to 7 eV. On the other 
hand, the PPP-CI results obtained using standard parameters, as also the earlier results 
obtained by Malloci et al .p predict that transition to the fifth excited state gives rise to 
the most intense peak. Thus, we conclude that, on the whole, the PPP-CI results calculated 
using the screened parameters are in better agreement with the experimental values than 


those computed using the standard parameters (cf. Table III). Screened parameter based 
calculations also predict that the wave functions of the excited states of the first five peaks 
(I-V) are dominated by single excitations (cf. Table 0. We also note that the TB model 
predicts the peak corresponding to the optical gap as the most intense one, located at 2.14 
eV, which is far away from the experimentally obtained value both in terms of peak location, 
and relative intensity. Therefore, we infer that the inclusion of electron correlation effects 
is essential for the correct quantitative description of the optical properties of graphene 
quantum dots. 
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Figure 5: Computed linear optical absorption spectrum for 16 atoms DQD, obtained by 
employing screened as well as standard parameters. In both the cases, the spectrum has 
been broadened with a uniform line-width of 0.1 eV. 

DQD-30 and DQD-^8 


In Figs. [6] and [7J we present the computed linear absorption spectra for DQD-30 and 
DQD-48, respectively. Information related to the energies, transition dipoles, and many- 
particle wave functions of excited states contributing to various absorption peaks for DQD-30 


are presented in Tables VI and VII of the Appendix 


As far as DQD-30 is concerned, our computed absorption spectrum has been compared 
with the experimental data of dibenzo[bc,kl]coronene (C 3 oH 14 ). 45 The experimental UV spec¬ 
trum obtained by Clar and Schmidt,® exhibits peaks at 2.55 eV and 3.61 eV. The position 
of the first peak at 2.55 eV is in good agreement with the computed value of optical gap 


at 2.08 eV (2.31 eV) (c/. Table III) obtained using screened (standard) parameters in the 


PPP-CI model. This peak is ^/-polarized, and corresponds to the li? 2 u state. The first ex¬ 
cited state of B 3u symmetry is dipole forbidden due to the particle-hole symmetry, and lies 
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Figure 6: Computed linear optical absorption spectrum for 30 atoms DQD, obtained by 
employing screened as well as standard parameters. In both the cases, the spectrum has 
been broadened with a uniform line-width of 0.1 eV. 


above the optical gap. Its energy is 2.25 eV (2.43 eV) obtained using screened (standard) 
parameters and is dominated by the \H — 2 —» L) + c.c. excitation. In addition, the ex¬ 
perimental peak at 3.61 eV agrees extremely well with the screened parameter x —polarized 
peak at 3.45 eV. This peak corresponds to a B^ u state, whose wave function is dominated 
by the \H — 2 —)• L) — c.c excitation. Thus, same excitations contribute to the wave func¬ 
tions of the first dipole-forbidden and dipole-allowed B^ u states, it is just that their relative 
signs are opposite due to the orthogonality constraint. Further, the most intense peak of 
the experimental spectrum is situated at around 6.20 eV, 61 which is in excellent agreement 
with the computed screened parameter value of the most intense peak at 6.28 eV. This peak 
also corresponds to a B% u state whose wave function consists of singly as well as doubly 


excited configurations (cf. Table VII). On the other hand, our standard parameter calcu¬ 
lations predict the most intense absorption peak at 5.10 eV corresponding to a state, 
whose wave function consists mainly of single excitations, dominated by the configuration 
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\H — 2 —> L + 2) (c/. Table VI). Thus, our computations imply that screened parameter val¬ 
ues are in better overall agreement with the experimental results than standard parameter 
values, and TB model predictions. While wave functions of the excited states corresponding 
to various peaks are dominated by single excitations, however, several states also derive 
significant contributions from the double excitations, hinting at the importance of electron 
correlation effects. 

In case of DQD-48, because of comparatively larger number of electrons in the system, 
the size of the MRSDCI calculations became excessively large. Therefore, we froze two 
lowest lying occupied orbitals, and deleted their particle-hole counterpart virtual orbitals 
which were highest in energy. With this approximation in place, the Cl problem reduced 
to that of 44 electrons, distributed over 22 occupied, and as many virtual orbitals, render¬ 
ing the calculation tractable. In order to benchmark this procedure, we also adopted the 
same methodology for DQD-30, and present the results of the calculations performed using 
screened parameters in Fig. [8] of the Appendix. It is observed that all the features of the 
optical spectra are preserved even after freezing the orbitals. However, the frozen spectrum 
is slightly blue-shifted as compared to the unfrozen one, with the corresponding changes 
being numerically acceptable. Next, we discuss our results for DQD-48 presented in Fig. [7J 
and Tables \VIU\ [TX} and |XJ 

We find that the first excited state of DQD-48 is a dipole forbidden B 3 U state, just as 
in the case of DQD-16, and is located at 1.38 eV (1.52 eV) as per our screened (standard) 
parameter calculations. Both the calculations predict it to be lower than the first dipole 
allowed state 1 B 2ui although the energy difference is much smaller as compared to the case 
of DQD-16. The wave function of this B% u state is dominated by the double excitation 
\H —> L + 1]H —> L) + c.c., and it will be of considerable interest if the future experiments 
on DQD-48 are able to locate this state relative to the optical gap. The first dipole-allowed 
peak in the absorption spectrum of DQD-48, corresponding to the li? 2 u state of the spectrum 
as in case of smaller DQDs, is computed at 1.61 eV (1.40 eV) based upon standard (screened) 
parameter based PPP-CI calculations. The wave function of this state is dominated by the 
| H —¥ L) excitation as in case of smaller dots, but it also derives significant contribution 
from the \H — 1 —y L + 1) configuration. The most intense peak of the absorption spectrum 
computed with the screened parameter is peak II corresponding to the 2B2 U state located 
at 2.19 eV, with the wave function dominated by the \H — 1 —y L + 1) configuration, 
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Figure 7: Computed linear absorption spectra (with a line-width of 
DQD, obtained by employing standard as well as screened 


0.1 eV) for 48 atoms 
parameters 


but also with a significant contribution from a triply excited configuration. On the other 
hand, standard parameter calculations predict peak Xll to be the most intense one, which 
is dne to a high-energy B- 2 U state at 5.54 eV, with the wave function dominated by several 
configurations. Such a large difference in the locations of the most intense peak predicted by 
standard and screened parameter calculations, can be easily tested in experiments. As far as 
general comparison between the standard and screened parameter calculations is concerned, 
besides the red shift of the screened parameter results compared to the standard ones, we 
find that only the first two peaks of the computed spectra have excited state wave functions 
which are qualitatively similar. We also note that compared to the two smaller DQDs 
discussed earlier, DQD-48 excited states exhibit significantly more contribution from doubly 
excited configurations. This trend implies higher contribution of electron-correlation effects 
in DQD-48. 
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IV. CONCLUSIONS 


In this work, very large scale correlated calculations employing the PPP Hamiltonian were 
carried out on diamond shaped graphene quantum dots of increasing sizes, namely, DQD- 
16, DQD-30, and DQD-48, and their optical as well as electronic properties were computed. 
Calculated linear optical absorption spectra of DQD-16 and DQD-30 were found to be in very 
good agreement with the experimental data of pyrene (CigHio) and dibenzo[bc,kl]coronene 
(C 30 H 14 ), which are their respective structural analogs with hydrogen passivated edges, thus 
justifying the essential correctness of our methodology. Some of the important conclusions 
we can draw from our correlated-electron calculations are: (i) the first peak corresponding to 
the optical gap is not the most intense, in contrast with the predictions of the tight-binding 
model, (ii) with the increasing size of the quantum dot, the absorption spectrum exhibits 
a red shift, (iii) the optical transition to the first excited state of B% u symmetry is dipole 
forbidden and it lies below the optical gap for DQD-16 and DQD-48, (iv) optical properties 
of the dots are sensitive to the projected corners of the system, therefore, they can be tuned 
by attaching suitable functional groups there. Thus, we hope that our work will spur further 
experimental activity in this field, so that our predictions on the excited states of DQD-48 can 
be tested in future experiments. Furthermore, recently Mullen and coworkers have stabilized 
graphene quantum dots with chlorine passivated edges) 58 Therefore, it will be of interest if 
chlorine passivated DQDs can be synthesized and their optical properties measured, so as to 
investigate the influence of the nature of edge passivation on the electro-optical properties 
of graphene nanostructures. 

In this work, we restricted ourselves to the study of linear optical properties of these 
quantum dots, but it will be quite interesting also to study the nonlinear optical response of 
these systems such as two-photon absorption, third harmonic generation, and photoinduced 
absorption. Calculations along those directions are underway in our group, and the results 
will be submitted for publication in future. 
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Appendix A: Influence of orbital freezing and deletion on the optical absorption 
spectrum of DQD-30 

In order to benchmark our orbital freezing and deletion approach aimed at reducing the 
size of the Cl calculations, in the figure below we present the results of optical absorption 
spectra of DQD-30 computed using all the electrons and orbitals, and by freezing (deleting) 
two lowest (highest) occupied (virtual) orbitals. Apart from a slight blue shift in the frozen 
orbital calculation, computed spectra are quite similar. 

Appendix B: Calculated Energies, Wave Functions, and Transition Dipole Moments 
Of The Excited States Contributing To The Linear Absorption Spectra 

Following tables represent the excitation energies, dominant many-body wave-functions, 
and transition dipole matrix elements of excited states with respect to the ground state 
(1 1 A 9 ). The coefficient of charge conjugate of a given singly excited configuration is ab¬ 
breviated as ’c.c.’, while the sign (+/-) preceding ’c.c.’ indicates that the two coefficients 
have (same/opposite) signs. For the doubly excited, and the higher order configurations, no 
dr/- sign precedes c.c. because more than one charge-conjugate counterparts are possible, 
each with its own sign. Label DF associated with a peak implies that the excited state in 
question is dipole forbidden. 
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Figure 8: Computed absorption spectrum (with a 
obtained by freezing orbitals (dotted line) and 


line-width of 0.1 eV) for 30 atoms DQD, 
without freezing orbitals (bold line). 
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Table V: Excited states giving rise to the peaks in the linear absorption spectrum of 16 
atoms DQD, computed employing the QCI {B 2u states) and MRSDCI ( B :iu states) 
approaches, along with the screened parameters in the PPP model Hamiltonian. 


Peak 

State 

E (eV) 

Transition Dipole (A) 

Dominant Configurations 

DF 

1 1 Bs u 

2.82 

0.000 

|H -> L + l)-c.c.(0.5730) 

\H —>■ L + 4; H — > L)+c.c.(0.1090) 

ly 

^ l B 2u 
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|H -> L) (0.8802) 

u x 

2 l B Zu 
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Illy 

3 1 B 2u 

5.37 
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77 — 1 -> L + 1) (0.8596) 

77 — 1 -> L + l;77-> L;H -> L)(0.1113) 

IVy 

&b 2u 

6.22 
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|77-2->L + 2)(0.6716) 

77 -> L + 5)-c.c. (0.3187) 

14 
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77 —L + 2; 77 — )• L)+c.c.(0.2700) 
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lO 1 ^ 
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VII X 
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n x B iu 
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|77 - 3 L + 4)—c.c.(0.2626) 


29 




Table VI: Excited states giving rise to the peaks in the linear absorption spectrum of 
DQD-30, computed employing the MRSDCI approach along with standard parameters in 

the PPP model Hamiltonian. 


Peak 

State 

E (eV) 

Transition Dipole (A) 

Dominant Configurations 

ly 

1 1 R2m 

2.31 

1.748 


|H -> L) (0.8223) 






77 — 1 —> L + 1)(0.1632) 

DF 

1 1 ^3m 

2.43 

0.000 


|77-2 -> L)-c.c.(0.4970) 





| H 

L + 1; 77 -> L)+c.c.(0.2365) 

HxSzy 


3.83 
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\H - 2 —> L)+c.c.(0.5432) 





| H 

— 5 —» 7,; 77 —» L)+c.c.(0.1563) 


4 1 S 2 „ 
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—^ 7/ + 1; 77 —>• L; 77 —)• L)(0.1507) 

Illy 
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- -1 
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1.047 

1 H 

- 1 -> L;77 L)—c.c.(0.3713) 
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| H 

-> L + 1; 77 -> L)-c.c.(0.2197) 
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1 5 1 B 2u 

5.70 
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| H - 

1 —)• L; 77 — 2 —> L)-c.c.(0.2148) 
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20 1 B 2u 

6.22 

1.113 


77 — 4 -> L + 4)(0.3458) 





| H - 1 

—>• L + 1; 77 —> L; 77 —> L) (0.2964) 
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Table VII: Excited states giving rise to the peaks in the linear absorption spectrum of 
DQD-30, computed employing the MRSDCI approach along with screened parameters in 

the PPP model Hamiltonian. 


Peak 

State 

E (eV) 

Transition Dipole (A) 

Dominant Configurations 

ly 

l 1 b 2u 

2.08 

2.036 

\H -» L}(0.8387) 

\H-l ->• L + 1)(0.1389) 

DF 

VB 3u 

2.25 

0.000 

\H-2-> L)+c.c.(0.4918) 

\H -»• L + 1; H L)-c.c.( 0.2691) 

II, 

VB 3u 

3.45 

1.443 

\H - 2 -»• L)-c.c.(0.5819) 

Illy 

4 1 B 2u 

3.71 

1.654 

\H- 1 -> L + 1)(0.7000) 

|ff-4-v i)+c.c.(0.2500) 

|ff - 1 -► L + 1; H -> L- H L)( 0.2466) 

n 4 

7 l B- iu 

4.23 

0.594 

\H - 6 ->• L)-c.c.(0.4901) 

|ff - 1 -> L; £7 -» L)+c.c.(0.2549) 

Vx&y 

11 1 b 3u 

4.90 

1.119 

|if-3->L + l)+c.c.(0.5083) 

|ff - 5 -► L- H -»• L)+c.c. (0.1442) 

|ff-7-»L + l;fT-» L)-c.c.(0.1350) 


10 1 b 2u 

4.83 

1.221 

|fi — 2 —>■ L + 2) (0.5634) 

|ff-2-»L + l;ff-» L)c.c.(0.2359) 

\H-1 ->• L + 5)+c.c.(0.2065) 

Vly 

15 1 B 2u 

5.32 

0.632 

if —>■ L + 3; if —>■ L)+c.c.(0.3245) 

ff — 7 —>■ L)—c.c.(0.3079) 

- 1 -» L; H - 2 L)—c.c.(0.2057) 

VII, 

20 1 B 3u 

5.69 

0.431 

|ff- 2 -> L + 4)-c.c.(0.4019) 

|ff — 1—>L+l;ff —>■ L + l)+c.c.(0.2542) 

|ff-l->L;ff-4-> L)+c.c.(0.2371) 

VIII, 

23 l B 3u 

5.97 

0.498 

|ff - 7 -> L + 2)+c.c.(0.3404) 

\H — >L + 2;ff — 3— >■ L)c.c.(0.3106) 

IX, 

29 1 B 3u 

6.28 

1.063 

if — 3 —> L + 5)+c.c.(0.3091) 

if — 6 —» L + 4)—c.c.(0.2390) 
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Table VIII: Excited states giving rise to the peaks I y to VII y in the linear absorption 
spectrum of DQD-48, computed employing the MRSDCI approach along with standard 

parameters in the PPP model Hamiltonian. 


Peak 

State 

E (eV) 

Transition Dipole (A) 

Dominant Configurations 

DF 

1 1 -P3« 

1.52 

0.000 

\H —> L + 1; H —> L)—c.c.(0.4873) 

|77-4->L)+c.c.(0.2313) 

Iy 

^ l B 2u 

1.61 

1.402 

|H ->L) (0.7467) 

|H - 1 -> L + 1)(0.3061) 

Ily 

3 1 i?2u 

2.57 

2.180 

|H - 1 -> L + 1)(0.5160) 

\H - 1 —> L + 1;77 -> L;77 -> L) (0.4413) 

77 — 2 —> L)+c.c.(0.2739) 

UI X 

3 1 i?3u 

2.76 

1.043 

77 -* L + 1; H —»• L)+c.c.(0.4903) 

77 -> L + 5)-c.c.(0.2202) 

|77^L + 1;77-1^L + l)+c.c.(0.1449) 

IVy 

5 1 B 2u 

2.88 

1.067 

77 — 2 —>■ L)+c.c.(0.4765) 

|77 - 1 -> L + 1) (0.2689) 

|T7 — 1 —> L + 1; 77 —> L; 77 -> L)(0.2549) 

V x 

B ?ya 

3.38 

0.933 

|77 — 4 ->L)-c.c.(0.4874) 

|77 - 6 -> L; 77 -> L)+c.c.(0.2237) 

|77 - 3 -> L + 1)—c.c.(0.1679) 

Vly 

§ l B 2u 

3.66 

1.727 

77 — 3 L; 77 —)• L) (0.3471) 

77 - 1 L + 1) (0.2908) 

VIIy 

10 1 b 2u 

4.16 

0.939 

77 — 2 —)• L + 2)(0.3609) 

|77 — 1 -> L + 7)+c.c.(0.3331) 
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Table IX: Excited states giving rise to the peaks VIII y to X II y in the linear absorption 
spectrum of DQD-48, computed employing the MRSDCI approach along with standard 

parameters in the PPP model Hamiltonian. 


Peak 

State 

E (eV) 

Transition Dipole (A) 

Dominant Configurations 

VUIy 

13 1 B 2u 

4.37 

0.596 

|T7 — 4—»L + 1;T7 —> 7)c.c.(0.3129) 

77 — 4 —> 7 + 4) (0.2685) 

|77 - 1 -> 7; H - 4 -> 7)-c.c.(0.2191) 

IXx&y 

1 s 1 b 3u 

4.59 

0.839 

\H - 3 -> L + 1)—c.c.(0.4035) 

\H — 8 —)• 7; 77 —> L + l)c.c.(0.1534) 

|H -> 7 + 5}-c.c.(0.1371) 


]& 1 B 2u 

4.65 

0.231 

|77-^7 + 4;77->7 + l)-c.c.(0.2687) 

\H - 6 -> L + l)+c.c.(0.2665) 

\H -> L + 9; H -> 7)+c.c.(0.2156) 

Xy 

20 1 R 2u 

4.80 

0.347 

\H — AL + 1;H 7)c.c.(0.2900) 

\H - 1 -» L + 5; H -> 7)c.c.(0.1942) 

77 — 2 —> 7 + 2) (0.1894) 

X IxSzy 

20 1 7? 3 „ 

4.96 

0.237 

\H — 9 -> L + l)+c.c.(0.4364) 

77 — 8—)• 7; 77 —> 7 + l)c.c.(0.1441) 

77 — 10 —> L; 77 —> L)—c.c. (0.1418) 


23 1 B 2u 

4.98 

0.634 

|77 — 2 ->7 + 2) (0.2775) 

|77 - 1 -> 7 + 5; 77 -> L)c.c.(0.2171) 

|77 — 4—>7 + l;77 —> 7)c.c.(0.2162) 

Xlly 

33 1 B 2u 

5.54 

1.183 

77 — 4 —>■ 7 + 4) (0.2954) 

|77 — 13 —> L)+c.c.(0.2185) 

77 — 1—»7 + l;77 —>■ 7 + 3)c.c.(0.2098) 


33 




Table X: Excited states giving rise to the peaks in the linear absorption spectrum of 
DQD-48, computed employing the MRSDCI approach along with screened parameters in 

the PPP model Hamiltonian. 


Peak 

State 

E (eV) 

Transition 

Dipole (A) 

Dominant Configurations 

DF 

1 1 R3u 

1.38 

0.000 

77 —> L + 1; H —> L)— c.c.(0.5106) 

|77 -> L + 3) - c.c.(0.2285) 

ly 

1 1 R2u 

1.40 

1.770 

|77 —>■ L) (0.7890) 

|77- l^L + l) (0.2338) 

Ily 

2 1 i?2u 

2.19 

2.910 

|77- l^L + l) (0.6593) 

|77 - 1 -> L + 1; H -> L; H -> L)(0.4735) 

U I, 


3.05 

0.984 

|77 — 3 -> L)+c.c.(0.4918) 

|77-1^L + 1;77^L + l)+c.c.(0.1791) 

n 4 


3.32 

1.135 

77 —>-L + l;77 — 2—> L)c.c.(0.3569) 

\H-1->L + 1\H->L + l)+c.c.(0.3481) 

Vy 

9 1 i?2« 

3.61 

1.313 

\H —> L + 4; H —)• L)+c.c.(0.3971) 

|77 - 9 -> L)+c.c.(0.2276) 

VI X 

16 1 B 3u 

3.96 

1.163 

77 — 4 —» L + l)+c.c.(0.3369) 

|i7->L + l;T7-2-> L)c.c.(0.2805) 

V Hy 

24 1 B 2u 

4.38 

0.661 

77—1—)-L + l;77—1—>L + 1;77—> L)(0.2485) 

77 — 2 —>■ L + 2)(0.2256) 

77 —> L + 1; 77 —»L + 1;T7 — 2—L)+c.c.(0.1943) 

V IH y 

33 1 B 2u 

4.78 

0.520 

77 — 1—)-L + l;77 —>L + 4)c.c.(0.3444) 

|77 - 1 -» L + 5; 77 -> L)c.c.(0.1999) 

|77 - 1 -> L + 10)—c.c.(0.1855) 
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